Transport properties of a two-dimensional electron liquid at high magnetic field 
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The chiral Luttinger liquid model for the edge dynamics of a two-dimensional electron gas in a 
strong magnetic field is derived from coarse-graining and a lowest Landau level projection procedure 
at arbitrary filling factors p < 1 - without reference to the quantum Hall effect. Based on this model, 
we develop a formalism to calculate the Landauer-Biittiker conductances in generic experimental 
set-ups including multiple leads and voltage probes. In the absence of tunneling between the edges 
the "ideal" Hall conductances {dj = e^u/h if lead j is immediately upstream of lead i, and Gij = 
otherwise) are recovered. Tunneling of quasiparticles of fractional charge e* between different edges 
is then included as an additional term in the Hamiltonian. In the limit of weak tunneling we 
obtain explicit expressions for the corrections to the ideal conductances. As an illustration of the 
formalism we compute the current- and temperature-dependent resistance Rxx{I,T) of a quantum 
point contact localized at the center of a gate-induced constriction in a quantum Hall bar. The 
exponent a in the low-current relation Rxx{I,0) ~ ja-2 gjjQ^g ^ nontrivial dependence on the 
strength of the inter-edge interaction, and its value changes as e*VH, where Vh ~ hl/ue^ is the Hall 
voltage, falls below a characteristic crossover energy hc/d, where c is the edge wave velocity and d 
is the length of the constriction. The consequences of this crossover are discussed vis-a-vis recent 
experiments in the weak tunneling regime. 



I. INTRODUCTION 

The Quantum Hall Effect (QHE) has been for the last 
twenty years an amazingly rich source of experimental 
a nd t heoretical resultsi^Si^ (for a review see also Refs. 
W^d). Exotic concepts such as incompressible quantum 
Hall liquids, fractionally charged quasiparticles^'^'^'^^'^^ 
and Composite FermionsiSiH have become part of every- 
day's language of physics. One of the most interesting 
developments triggered by the QHE has been the real- 
ization that the edge of a quantum Hall liquid^'* pro- 
vides a clean realization of the chiral Luttinger Liquid 

As is well known, the Luttinger liquid (LL) concept 
- first introduced by Haldane^^, building on an earlier 
exact solution of the Luttinger liquid modeliSiiS - is the 
accepted paradigm for the low-energy behavior of inter- 
acting Fermi liquids in one dimension. In the LL model 
two types of fermions - right movers and left movers - 
are coupled by an interaction of strength g. Each type 
of fermion by itself forms a chiral Fermi liquid, and its 
density fluctuation Spa{x) {a = left (L) or right (R)) can 
be expressed as the derivative of a bosonic displacement 
field 4'a{x) that satisfies the commutation relations 

[(paix), 'f)f3{x')] = iTTSaSaf3Sig-n{x - x') , (1) 

where Sa — 1 for a = L and Sa — —1 for a — R. The in- 
teraction between right and left movers can be eliminated 
by a transformation (canonical up to a scale factor) that 
preserves the relation between the net current and the 



displacement fields and leads to two independent chiral 
fields (/)^ which now satisfy the anomalous commutation 
relation 

[4>'aix), (f>'p{x')] = iTTe~'^^SaSa/3signix ~ x') , (2) 

where & — tanh^^ | is a measure of the strength of 
the original left-right coupling. The xLL model arises 
when one considers just one of these two fields, with the 
commutator (0). 

The anomalous commutator leads to a rich phe- 
nomenology, including absence of the usual electron 
quasiparticles, anomalously slow decay of correlation 
functions, nonlinear transport properties etc... Needless 
to say these effects are very difficult to observe experi- 
mentally, due to the dramatic impact of even a modest 
concentration of impurities on the properties of a one- 
dimensional quantum system^^i^^i^^i^^i^^i^^i^^i^'''. 

It was therefore welcome news when, in a seminal 1990 
paper, WeriiSiiS*S& showed that the density fluctuation 
excitations at the edge of an incompressible quantum 
Hall liquid at filling factor v = l/q {q—odd integer) cor- 
respond to those of a xLL with e~^^ = v. Unlike one- 
dimensional metallic systems, the edge of a quantum Hall 
liquid is essentially unaffected by disorder, so the xLL 
ideas could finally be put to an accurate experimental 
test^Si^. Following Wen's insight the analysis was ex- 
tended to more complex hierarchical QHE states, where 
it turned out that one can have multiple branches of edge 
excitations (i.e., multiple xIAjs), some propagating in op- 
posite directions, and disorder plays a role in ensuring the 
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correct value of the quantized Hall conduct ance^ii^. 

Understandably, these papers created a widespread be- 
lief that the xLL behavior of the edge is inextricably tied 
to the QHE in the bulk. For one thing, the energy gap of 
the quantum Hall liquid state was believed to be essential 
to ensure that the low energy excitations are confined to 
the edges of the system. It thus came as a big surprise 
when Grayson et al.^ reported that the xLL could be 
observed in a whole range of filling factors j < ly < I 
and was apparently unrelated to the quantization of the 
Hall conductance. 

In Section II of this paper we will argue that the valid- 
ity of the xLL model for the edge dynamics of a two- 
dimensional electron liquid at high magnetic field fol- 
lows from elementary semiclassical considerations, which 
should be valid at any filling factor and have nothing to 
do with the occurrence of the quantum Hall effect. The 
essential point is that the hydrodynamic modes, obtained 
by "integrating out" fluctuations that are rapidly vary- 
ing in time (on the scale of the cyclotron frequency) and 
in space (on the scale of the magnetic length)^"*, are au- 
tomatically bound to the regions of space in which the 
gradient of the equilibrium density differs from zero, i.e. 
to the edges of the system. In addition, the algebra of the 
edge density fluctuations (precisely defined in the next 
section) follows from the algebra of the projected density 
operators, when the latter is averaged on a length scale 
that is large compared to the magnetic length. 

In Sections III-IV we develop the formalism for the cal- 
culation of the Landauer-Biittiker (LB) conductances^ 
for generic experimental arrangements including multi- 
ple terminals connected to the system by leads. Ordinar- 
ily, the LB theory expresses the conductance Gij (which 
connects the current in the i-th terminal to the voltage 
applied to the j-th one) in terms of the transmission prob- 
ability of an electron quasiparticle from one terminal to 
the other. But, in the present case, there are no electron 
quasiparticles. Instead, the voltage applied to one termi- 
nal induces a train of collective waves which propagate 
along the edges of the system and eventually feeds a cur- 
rent into several different terminals. It is not surprising 
therefore that the conductance can be expressed solely 
in terms of the displacement field propagators along the 
edges of the system. We show that this approach, in the 
absence of inter-edge coupling, yields the ideal Hall con- 
ductances even in the presence of inhomogeneities that 
cause partial reflection of edge waves. 

Deviations from the ideal Hall effect can and do oc- 
cur when the possibility of inter-edge tunneling is taken 
into account. This subject is taken up in Sections V- 
VI. Due to its quantum mechanical origin tunneling is 
not included in the semiclassical hydrodynamic descrip- 
tion and must be introduced "by hand" . We describe 
tunneling in terms of two parameters, the tunneling am- 
plitude r and, most importantly, the charge e* of the 
quasiparticles that are transferred from one edge to the 
other. Since a fundamental theory of e* at general fill- 
ing factors is not yet available one may choose to treat 



e* as a phenomenological input parameter, whose value 
may be determined from experiments. Alternatively one 
can choose e* — ve which is believed to be correct at 
V — 1/q, where q is an odd integer. In terms of F and 
e* we can finally calculate the corrections to the ideal 
Hall conductances: the final expressions involve the dif- 
ferential tunneling conductance, i.e. the derivative of the 
tunneling current with respect to the potential difference 
between the two edges. 

In Section VII we present a perturbative study of the 
nonlinear resistance RxxiI,T) of a quantum point con- 
tact situated within a constriction in a quantum Hall 
bar—-. The perturbation theory is valid for R^-x <C 
This study generalizes Wen's original treatment of this 
phenomenonSi and the later study by Moon and Girvin^S 
by including the effect of an inhomogeneous short-ranged 
inter-edge interaction, i.e., an interaction that is strong 
in the region of the quantum point contact, but becomes 
weak as one moves away from it. In Wen's paper a repul- 
sive, but translationally invariant, inter-edge interaction 
leads to a decrease in the tunneling exponent a defined by 
Rxx{I, 0) ~ or Rxx{0, T) ~ T"-2. This would make 
the behavior of the resistance even more singular than in 
the theory without inter-edge coupling at low tempera- 
ture and bias voltage. Our calculations indicate that the 
interplay of the inter-edge interaction with the broken 
translational invariance alters the relationship between 
the tunneling exponent and the strength of the interac- 
tion in the constriction region. The new relationship is 
relevant when either e*l/f/, {Vh = ^ being the Hall 
voltage) or fcsT are above a geometric energy scale ^, 
where c is the edge wave velocity and d is the length of 
the constriction. For realistic values of the parameters 
this energy scale is in the range of 100 mK. Above this 
"crossover" energy the tunneling exponent turns out to 
be larger than expected from the noninteracting theory 
and a fortiori, from Wen's interacting theory. 

All these results suggest that a quantitative compari- 
son between theory and experiment cannot ignore the in- 
teractions between the edges of the quantum Hall liquid 
in the region of the constriction. In particular the ex- 
ponents of the current-voltage relationship may be non- 
universal in the experimentally accessible range of tem- 
peratures, reverting to universal values only at extremely 
low temperatures. Evidence for nonuniversal behavior in 
the tunneling exponents has recently surfaced from sev- 
eral different points of vie vj^^i'^°i'*-^i'^^i"^^i'^'^ . 



II. DERIVATION OF THE CHIRAL 
LUTTINGER LIQUID MODEL 

Consider a two-dimensional electron liquid in a strong 
perpendicular magnetic field B — —Bz such that all the 
electrons reside in the lowest Landau level (LLL). The 
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hamiltonian, projected within the LLL, has the form 
H = \ dr dr'p{r)V{r-r')p{r') 



drVo{r)p{r) 



(3) 



where p{r) is the number density operator projected in 
the LLL, V{r — r') is the electron-electron interaction 
potential, and Vb(r) is an external potential. Both the 
kinetic energy and a self-interaction-removing term are 
just constants, and have therefore been dropped. 

Next we write the density operator as the sum of 
the classical equilibrium density po{v) and a fluctuation 
Sp{ry. 



p{r) = pa{r)+Sp{r) , 
where Po{y) is determined by the equation 



dr'V{r~r')po{r') + Voir)= p 



(4) 



(5) 



and /i is a constant fixing the total number of particles. 
This gives (again, up to a constant) 



dr / dr'Sp{r)V{r-r')Sp{r') . 



(6) 



The commutation relations between projected density 
operators at different r are easily deduced from the well 
known result^ 



[p(q),p(k)] = (e'=* 



-kq'l^/2 



p(k + q) 



(7) 



1/2 

where k = k^ + iky, q — Qx + iqy, (■ = {j§) is the mag- 
netic length, and {x, y, z) form a right-handed coordinate 
frame. 

Since we are interested in the dynamics of long wave- 
length density fluctuations, we expand the right-hand 
side of to leading order in k£, q£ and transform to 
real space. This gives 



[p(r),p(r')] ~ ife,jdiPo{r)djS{r 



(8) 



where i,j denote cartesian components in the {x,y) 
plane, di = Cij is the two-dimensional Levi-Civita 
tensor, and repeated indices are summed over. Notice 
that we have replaced the density operator p{r) on the 
right hand side of Eq. (jSj) by its equilibrium expectation 
value po(r): this is legitimate as long as we are interested 
only in the linear dynamics of small fluctuations about 
the equilibrium state. 

The remarkable feature of Eq. © is that the commu- 
tator is proportional to the derivative of the ground state 
density. This implies that hydrodynamic density fluctua- 
tions are bound to regions where the equilibrium density 
varies, and are absent from the regions of constant den- 
sity. This can be seen most clearly by writing down the 



equation of motion for the density fluctuation, which is 
easily seen to have the form 

dt5p{v,t)^f{d,po{v))e,,djj dv'V{v,v')6p{v',t) . (9) 

Because this equation agrees with what one finds by tak- 
ing the large magnetic field limit of the hydrodynamic 
Euler equations^ we will call our approach "hydrody- 
namical" . 
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FIG. 1; a) Description of the device used in the experiments. 
A two-dimensional electron gas is in contact with several 
reservoirs. Two, the source (S) and the drain (D), are used 
to inject and extract a current. The others (four in the fig- 
ure) are used as voltage probes. Each reservoir contacts two 
edges (labelled by roman numerals). Gate voltages (G) cre- 
ate a depletion zone and force the edges to stay close. The 
current flows are shown in the figure, b) Expanded view of 
the edge region. Density fluctuations exist only in a limited 
region around the x axis (the shaded region in the figure). 
We integrate over the direction y perpendicular to the edge 
to obtain an "edge density fiuctuation" that depends only on 
the position along the edge. 

The fact that the hydrodynamic density fluctuations 
are proportional to the derivative of the equilibrium den- 
sity implies that they are concentrated near the edges of 
the system where the density profile has a strong varia- 
tion. Let us consider, for definiteness, the model depicted 
in Fig. ^ where the boundary of the electron liquid is 
divided by leads into different "edges" labelled by roman 
numerals. The leads are connected to reservoirs (labelled 
1-4 in Fig. 0). There are also two special reservoirs, the 
"source" (S) and the "drain" (D). Each reservoir is con- 
nected to two edges. We introduce a one dimensional 
coordinate Xa that keeps track of the position along edge 
a {a = i,ii,iii... in the example) growing continuously 
along the direction of the arrows, i.e. from source to 
drain. We denote by Xai the point where the lead com- 
ing from reservoir i contacts the edge a. At each point 
along an edge we attach a local y-axis normal to the edge. 
The density varies rapidly as a function of y and slowly 
as a function of x. We therefore introduce an integrated 
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edge density fluctuation 

Sp{xa) ^ I dy Sp{xa,y) 



(10) 



where the edge is located at about y ~ and the integral 
over y extends far enough to include the whole region 
in which the density fluctuations differ from zero (see 
Fig. From the algebra ||SJ) one can then derive the 
commutation rules for edge density fluctuations. In the 
Appendix El we show that 



■ l^{^a)Sa roc-/ 

^ Oal3(Jx^(>(Xo 

ZTT 



(11) 



where Sq = 1 if a is a left edge, Sq, = — 1 if a is a right 
edge, ly'lxa) is the equilibrium fllling factor (deflned as 
^{x) = 2ttPpo{x)) in the bulk contiguous to the edge 
labelled "a" . Edge fluctuations on different edges com- 
mute. 

For constant f Eq. Hll|) reproduces the Kac-Moody 
current algebra for the density fields in the chiral Lut- 
tinger liquid model. Due to the presence of the nonin- 
teger factor ly these commutation rules imply that each 
edge with v < 1 exhibits xLL behavior, even if all in- 
teractions between different edges are turned off. Notice 
that our derivation of the xLL has nothing to do with 
the quantum Hall effect: it only depends on the coarse- 
graining (the hydrodynamic approximation) and on the 
high magnetic field limit (the projection in the LLL). The 
derivation is valid for arbitrary value of v, whereas the 
quantum Hall effect occurs only at special values of v for 
which, as will be argued below, the strength of tunneling 
between different edge states becomes negligible. 

In terms of edge density fluctuations, the hamiltonian 
becomes 



H = 



Q/3 



dXn 



dx'i3Sp{Xa)V{Xa,x'i3)Sp{x'p) (12) 



where V{xa,x'p){= V{x'i^,Xa)) is constructed from the 
original interaction V{r,r') by putting r at position Xa 
along the edge a, and r' at position x'^ along the edge f3. 
By using Eq. (|ll|l and the hamiltonian H12|) , the equation 
of motion for edge density fluctuations is immediately 
found to be 



dtSp{xa) 



2^ 



dx'f3 dx^V{xa,x'p)6p{x'i3) 



(13) 

Rather than pursuing the solution of Eq. I|13|) in general, 
we shall henceforth restrict our attention to the special 
case in which all the edges share the same bulk density, 
i.e. i^ixa) = 1^ independent of a and x. It is convenient 
to define the "displacement field" (l){xa,t) such that 



6p{xa,t) = dx4>{xa,t). 
These fields satisfy the commutation relations 

[4'{Xa),(t>{x'f^)] = i^SaSign(x - x')6al3 



Assuming a time-dependence of the form 

^ix^,t) = 4,ix^)e-''^\ (16) 
we see that Eq. (|13|) takes the form 



iujdx^(j){xa) = / dx'i3dx^V{Xa,x'ij)d.j;'^4>{x'i^) 



The associated eigenvalue problem 



(17) 



iLodx^ifixa) ^ I dx'pdx^V{xa,x'p)dx>^ip{x'p) 



is hermitian, and has the following properties: 
1. All the eigenfrequencies LOn are real. 



(18) 



2. If Lpn{xa) is an eigenfunction with frequency uin 
then ip'l^{xa) is an eigenfunction with frequency 

3. The eigenf unctions (pna{x) form a complete basis 
in the Hilbert space with the completeness 

-i^sign(cJ„)93„(Xa)(p*(a;^)(9a;j^ = Sa5apS{Xa - x'^). 
n 

(19) 

and the orthonormality conditions 

^«Sq / dXa fn{Xa)dx^fm{Xa) ^ Slgn{uJn)5nrn (20) 
a •' 

The proof of these relations is provided in the Appendix 

m 

It is straightforward, with the help of the above rela- 
tions, to show that the density fiuctuation field can be 
expanded on the basis provided by the </9„'s as follows 



' ri>0 

(21) 

where n > specifies that only the positive frequency 
eigenf unctions are included in the sum and the &„s - one 
for each n > - are boson operators obeying the stan- 
dard commutation relation [&n,^|i'] = ^nn' ■ At the same 
time, the hamiltonian (|12|l takes the form 



(22) 



ri>0 



It is instructive at this point to solve the eigenvalue 
equation H18|) in a simple case. We consider just two 
parallel edges in a translationally invariant Hall bar ge- 
ometry. Of = 1 for the left edge and a = 2 for the right 
(14) edge (see Fig. The interaction is assumed to have the 
form 

m ">.-4)^(;;;j5:5;!ll!S:g) p« 
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where Vi and V2 are translationally invariant interactions 
between density fluctuations on the same edge and on 
different edges respectively (The coordinates xi, X2 are 
set up to have the same value on points at the same 
"height" on the two edges). We seek the solutions of 
Eq. in the form (p{xa) — ipa{k)e^^^'' . This leads to 
the 2x2 eigenvalue problem 







(P2 



ikv 
27 



Vl{k) V2{k) 

V2{k) Vi{k) 



^2 



(24) 

where Vi (fc) and V2 (k) are the Fourier transforms of Vi (x) 
and V2{x) (the upper part of the spinor refers to the left 
edge). The eigenvalues are 



LOk^±^^V,Hk)~Vi{k)\k\ , (25) 

and for each positive frequency there are two solutions: 
the "up-moving" one is 



1 



/kL 



UkC 

-Vke"' 



ikxi 



with fc > 0, and the "down- moving" one is 



1 



-UkC 



-ikxi 
-ikx2 



(26) 



(27) 



also with fc > 0. Here, as usual, we have normalized 
the eigenf unctions with the factor 1/\/L where L is the 
length of the edge. This length is assumed to be arbi- 
trarily large, and will not enter the physical results. On 
the other hand the presence of the normalization fac- 
tor 1/Vk is imposed by the orthonormality conditions. 
Since these conditions require u'j. 
Uk = cosh6'fe, Vk = sinh6'fe and 



— 1, one can write 



tanh2(?fe = 



V2(k) 
Vi{k) 



(28) 



Thus, we have recovered the standard expressions for the 
dispersion of the edge waves in the ordinary (nonchiral) 
Luttinger liquid model. However, we emphasize that, in 
the present model, the xLL behavior persists even if the 
interaction V2 is turned off. This is due to the anoma- 
lous commutator between density fluctuations on the 
same edge. 



J 



III. FORMULATION OF TRANSPORT 

In transport theory we need to calculate the change 
in the current li that flows into reservoir "i" due to a 
change in the potential Vj of reservoir "j" : 



(29) 



(The conductance matrix elements will in general de- 
pend on the initial values Vi of the applied voltages: these 
initial values will be referred to as "bias voltages".) The 
current will be considered positive when it enters a reser- 
voir and negative when it leaves it. Due to gauge invari- 
ance and current conservation the G^s (under steady- 
state conditions) satisfy the constraints 



(30) 



These constraints specify the values of the diagonal con- 
ductances Gii once the off-diagonal ones with i ^ j are 
known. 

The form of the edge current-density operator I{xa) is 
dictated by the continuity equation 



edtSp{xa) = dxl[xa) 
which immediately gives 

Ia{x) = edt(f>aix). 



(31) 



(32) 



This current density is positive when it flows along the 
direction of the arrows in Fig. ^ negative otherwise. 
Thus, the current flowing into terminal i is given by 



SI, 



(33) 



where 



- 1 if a enters i 
-1 if a exits i 
otherwise. 



(34) 



The linear response of the current density to a periodic variation of the electrical potential SV^xp) is given by 
Slix^) = '-f-T.J S^i^'p) dt([9*0(x„,t),9,;,0(x;3)])e'("+"')* , 



(35) 



where to is the frequency and (...) denotes the equilibrium average. A first integration by parts with respect to time 
gives 





51{xa) ^ i^z2 '^^'p ['^(^x'i^^{xa,x'i3;uj)]SV{x'p) , 
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where 



is the retarded displacement-field propagator, whose explicit expression in terms of "phonon eigenfunctions" is 



D{xa,x'p;uj) = I'X! 



ri>0 



Ll! — LOn + if] 



IT] 



(37) 



(38) 



In doing the integral by parts we have exploited the fact that dx' D{xa,x'p]t = 0+) = —2'Ki[(j){xa)Tdx'^<f>{x'p)\ oc 
Saf3S{xa — x'p) vanishcs unless Xa and x'^ coincide. It will be shown below that this condition is always satisfied in 
the relevant region of integration. 



A potential change SVj applied to the j-th lead can be 
modelled as a change of the potential on the two edges 
that enter and exit the reservoir. The change in poten- 
tial is considered uniform over the portions of the edges 
that run inside the leads, and drops to zero at the points 
of contact between the leads and the system (It must be 
borne in mind that what we are modelling here is the 
externally applied potential, not the full screened poten- 
tial that will appear all over the system in response to 
the external perturbation). Thus we see that the poten- 
tial change associated with reservoir j is described by the 
equation 

dxpSV{xp) = ^£,f3jS{x0 - x,3j)SVj , (39) 
j 

where the "contact functions" ^/3j are defined in Eq. (|34|l . 

We now combine Eqs. and ((23 ■ The integral 

over x'p can be immediately carried out (by parts) under 
the reasonable assumption that the phonon eigenfunc- 
tions decay exponentially for x —^ ±oo i.e., well inside 
the reservoirs. This is physically expected to happen as 
the one-dimensional edge channels broaden into a three 
dimensional reservoir. Mathematically, one must make 
sure that the eigenfunctions used to calculate the dis- 
placement propagator satisfy this boundary condition. 
The final result for the current arriving at reservoir i via 
edge channel a due to a potential disturbance applied to 
edge channel P by reservoir j with i j is 

( \ 

(40) 

The quantity within the round brackets is, by defini- 
tion, Gij. Note that this equation specifies only the 
off-diagonal elements {i ^ j) of the conductance matrix. 
This guarantees that Xai is macroscopically distinct from 



xpj and validates our integration by parts with respect 
to time. Diagonal elements Ga are determined by the 
continuity conditions 1)30(1 . 

As a simple example consider the calculation of the 
propagator in the translationally invariant geometry with 
short range interactions, so that the phonon eigenfunc- 
tions are labelled by a wave vector k and Uk = ck where 
c is the velocity of the edge mode. To ensure that the 
phonon eigenfunctions vanish for — s- oo we shift the 
wave vector k infinitesimally into the complex plane, set- 
ting k — > k + ir]Sgn{x) in Eq. ((26(1 and k —> k — ir]Sgn{x) in 
Eq. 1(27(1 . Next, we substitute these eigenfunctions into 





D 




' 


X X ^ 

-*1.3 -*!I,2 

I II 

^I.l ^11,1 


r 






1 1 

s 





FIG. 2: Schematics of a translationally invariant two-terminal 
device. 



Eq. ((38(1 and convert the sum over n into an integral along 
the real axis of the complex variable k. We readily find 
that only the poles at k — ±{uj + ir])/c contribute to the 
integral which thus yields 



i Ihn UjD{XajX'jj;Uj) — vQ^Xa — x'a) 



-uv 

,,2 



+ V<d{x'g 



-uv 



(41) 
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where u, v are the fc — > Hmits of Uk and Vk. 

Notice that in the case of decoupled edges {u = 1, 
V — 0) one has only upward propagation on the left edge 
and downward propagation on the right one. This makes 
the conductance Gij vanish unless the reservoirs j is "up- 
stream" of reservoir i, consistent with the definition of 
an ideal quantum Hall system^. It is straightforward, 
at this point, to compute the two-terminal conductances 
Gi2, G21 of the simple device shown in Fig. [3 Since the 
source and the drain reservoirs contact both edges, and 
^Qi = — for each edge, Eq. l|in|l gives us 

G12 = 4^e-2« ^ G21 . (42) 
h 

Interestingly, the presence of the factor e~^® — {u — t;)^ 
in the relation between the current and the source-drain 
potential does not imply a deviation from the ideal Hall 
conductance, since the relation between the Hall voltage 
(as measured by ideal voltage probes applied to the two 
sides of the Hall bar) and the source-drain potential is 
also modified by the same factor—. 

IV. REFLECTION AND TRANSMISSION OF 
EDGE WAVES 

Before proceeding to the calculation of the conduc- 
tances in the presence of inter-edge tunneling we wish 
to take a closer look at free edge waves in the presence 
of a constriction that breaks translational symmetry (see 
Fig. 13). This constriction can be created by depleting a 
portion of the sample by applying a voltage to metallic 
gates on top of the mesa. When an edge wave of finite 
wave vector k impinges on the constriction it is partially 
reflected and partially transmitted. How this affects the 
conductance depends crucially on the behavior of the re- 
flection coefficient r{k) in the limit A; ^ 0. If r(0) = 
then there is no correction to the ideal conductance (in 
the absence of tunneling); otherwise there will be one. 

To keep the analysis simple, we now assume that both 
Vi and V2 are short-ranged on the scale of the density 
variations: this means, in particular, that only points at 
the same value of x on opposite edges interact and Vap (x) 
has the form 

where 

r V2.1 X < -d/2 

V^ix) = { V2.2 \x\ < d/2 . (44) 

[ 1/2,3 X > d/2 

Since the potential has a step-like behavior, with three 
different values in the three regions 1, 2, and 3, we seek 



r 
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FIG. 3: The two types of constrictions we consider. On the 
right, a semi-infinite constriction. On the left, a more realistic 
constriction localized in a finite region of the sample. For 
clarity the lateral voltage probes are not shown. 

the solution in a piece- wise form. As in the standard scat- 
tering theory we label the full solution with the quantum 
number of the incident wave. For instance for an incident 
wave from the bottom on the left edge (^^^ (x) we seek the 
"up-moving" solution in the form 

r cpl{x)+r"ipi^{x) x<-d/2 
(x) = <^ A"^l (x) + (x) \x\ < d/2 (45) 

[t-iplix) x>d/2 

The wave vectors fci, ^2, and ks in regions 1, 2, and 
3 respectively are determined by the condition that the 
energy of the wave is not changed in the scattering, i.e. 

cifci = C2fc2 = csks (46) 

where ci, C2, and C3 are the sound velocities in the corre- 
sponding regions. We remark that the wave function tp^^ 
is labelled with the wave vector ki of the incident wave 
it originates from. In similar way one can construct the 
"down-moving" solution, 

(t^^tM x<^d/2 
^i{^) = { A%l^(x) + B'^iplix) \x\<d/2 , (47) 
y^i^{x)+r''^l^{x) x>d/2 

Notice that the spinor-like eigenfunctions Lp^'^\x) (see 
Eqs. (|26|) and (|27l) ') are those appropriate for each region. 
^u(d) ^ ^u(d) jQ^^ reflection and transmission amplitudes 
for the up- (down-) moving wave. The matching condi- 
tions are dictated by the physical requirement that there 
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is no accumulation of energy at the interfaces. This is the coefficients A, B, t and r can be determined. The 
equivalent to the requirement of continuity of the solu- solution, expressed in terms of the mixing angles, is 
tion at a; = ±d/2 and gives four conditions from which 



cos(/c2C?) cosh(6'i — 63) — i sm{k2d) cosh(2(?2 — &i — ^3) V C3 ' 
cos(fc2(i) sinh(0i - ^3) + i sin(fc2d) sinh(202 - 9i - O3) 
cos(fc2rf) cosh(0i — ^3) — i sin(fc2(i) cosh(2(?2 — Oi — 63) ' 



cosh(03 — 02)e 



(48) 



cos(fc2d) cosh(0i — 63) — i sin(fc2d) cosh(26'2 — 9i — 63) y C2 ' 
sinh(6'3 - 6*2)6'^^^'^ fE^ 



cos(fc2(i) cosh(6'i — 6*3) — i sin(fc2(i) cosh(2(?2 — 9i ~ O3) \ C2 



The expression for the coefficients in the "down-moving" 
solution can be obtained with the substitutions 



B"; B" 



and 



6*1, ci 



C3- 



A"^ (49) 



(50) 



The reflection and transmission coefficients satisfy the 
relation |rp + (ci/c3)|tp = 1 which follows from the con- 
servation law derived in Appendix El By sending d to 
zero one can examine the case described in the left panel 
of Fig. El When d = 0, one gets^ 



1 



cosh(6ii - 63) 
sinh(6li - 63) 



(51) 



cosh(0i 



3)' 



In Ea. (|42(l we have seen that the inter-edge interaction 
renormalizes the two terminal conductance G12 with the 
factor e~^^ which gives rise to an effective filling fac- 
tor V — ve^^^ . From this point of view by introduc- 
ing the effective filling factor in the different regions 
Vi — ve~^^' [i = 1, 2, 3), one can rewrite (|51|l aa^ 



t : 



2i>3 jci 



C3 



(52) 



We can now ask how the reflection of edge waves mod- 
ifies the conductances obtained for the translationally in- 
variant case at the end of the previous section. To keep 
the discussion as simple as possible consider first the situ- 
ation in which the inter-edge interaction is present only in 
the constriction region —d/2 < x < d/2. Four reservoirs 
are attached to the system above and below the constric- 
tion. A straightforward calculation with eigenfunctions 



of the form (|45|l and H47|l gives the following expressions 



for the dimensionless conductances 



hG,j . 



521 



1 



27ri 
931 = 

1 



941 



2m 



ik{x2-xi) 

dk— — -r(fc) 

, k-uj/c-iQ+ ^ ' 



-ik{x4,+xi) 

dk- r»-"(fc) 

, fc-w/c-iO+ ^ ' 



(53) 



with the transmission and reflection coefficients given by 
Eq.(|3HI) evaluated for 0i ^ 63 ^ 0. They read 



r(fc) = 



icos(cifcd/c2) -I- sin(ciA:d/c2) cosh(202) 



r"(fc) = - 



sin(cifcc?/c2) sinh(202) 



(54) 



s{cikd/c2) + sm{cikd/c2) cosh(202) 



The key observation at this point is that the exponential 
factors in Ea. (|53|) force to close the integration contours 
in the upper half-plane. Apart from the pole at k — 
Lu/c + iO+ there are no other poles in the upper half- 
plane since both t"(fc) and r"(fc) have poles in the lower 
half-plane. Thus, we obtain 



G21 - ^t^icj/c) 
h 

G41 - 4^r"(c./c) 
h 



(55) 



and in the limit of zero frequency one recovers the ex- 
act quantization that characterizes the ideal fractional 
QHE. This result can be understood by observing that in 
the long wavelength limit the constriction becomes fully 
transparent to the current. The situation is quite differ- 
ent in the case of the semi-infinite constriction where r(0) 
acquires a finite value: in this case one finds deviations 
from the ideal Hall conductance. It is amusing to see 
that the expressions (|55|) are similar in form to what the 
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Landauer-Biittiker theory would predict for a situation 
in which particles are physically backscattered from one 
edge to the other with probability r(0). However, up to 
this point, our theory does not allow for the transfer of 
charge between the edges. 



V. THE TUNNELING HAMILTONIAN 

It is now time to consider the effect of charge tunneling 
between different edges. The physical origin of tunnel- 
ing lies in the fact that the electron quasiparticles are 
not 100% localized on one or the other edge: the density 
matrix p(r,r') = (r)*iE'(r')) ('l'^(r) is the creation op- 
erator of a quasiparticle at position r) has a finite value 
even when r and r' are on different edges. This is true at 
all filling factors but, of course, the range of the density 
matrix depends dramatically on whether there are ex- 
tended quasiparticle states at the chemical potential, and 
this in turn depends on the filling factor. The fractional 
QHE is believed to arise at electronic densities such that 
there are no extended quasiparticle states at the chemi- 
cal potential, so that /o(r, r') is exponentially small when 
r and r' (on different edges) are separated by a distance 
larger than ^ I. This means that there is essentially no 
tunneling between the edges. Even in this case, however, 
tunneling can be induced by pushing two edges together 
as in the constrictions studied in the previous section. 

Since the physics of tunneling is lost in the hydrody- 
namic approximation, (which is local in space and there- 
fore does not allow for any direct connection between the 
edges) we need to put it back in the hamiltonian "by 
hand" . To this end we define a quasiparticle operator 
^J,(xa) which adds a charge e* (not necessarily equal to 
the electron charge e) localized at position Xa along the 
a edge. This is accomplished by requiring that '^]^{xa) 
satisiy the commutation relation 



-5ap5{xa- xL)^!l{xa). (56) 



We hasten to say that we do not know, in general, what 
the correct value of e* is. In some special cases, for exam- 
ple at filling factors of the form v = 2n+i with integer n, 
it is widely believed that e* = z/e, but there is no general 
theory for arbitrary filling factors. Let us then treat e* 
as a phenomenological parameter, and note that Eq. H56|) 
is satisfied bjiS 



^aiixc) =f7lexp 



X exp 



n>0 

-i—\ —Sa (p,iiXa)bn 
n>Q 



(57) 



where is an operator that commutes with all the 6's 
and 6^'s and increases the total charge Qa on the edge 
by-e* 

[UlQ0]=e*d^pUl. (58) 



The statistics of the quasiparticle is determined by the 
charge. If e* — e and v = 1 the fermion commuta- 
tion relations are satisfied, otherwise the quasiparticle 
has a fractional charge and a fractional statistics. An- 
other point to be made is that the creation of a fractional 
charge is not in contradiction with the quantization of 
the electric charge: the fact that the charge on the edge 
varies by a fractional amount simply means that a com- 
pensating fractional variation must be occurring deep in 
the reservoirs to keep the total charge in the universe an 
integer. 

In terms of the quasiparticle operators, the tunneling 
between two edges ( "left" and "right" ) coupled by a con- 
striction at a; = is described by the tunneling hamilto- 
nian 

Ht = T: *l(0)*fl(0) : +r* : ^'Jj(0)*l(0) :, (59) 

where F is a (phenomenological) tunneling amplitude and 
: . . . : indicates the normal ordering. One can of course 
consider more general situations in which tunneling oc- 
curs simultaneously at different points. 



VI. THE TUNNELING CONDUCTANCE 

Let us consider the case of just two edges coupled by 
a constriction at x = 0. The complete hamiltonian is 



(60) 



n>0 



where the tunneling term Hti given by Eq. (|59|l . intro- 
duces an interaction between the formerly free bosons 6„. 
At the same time, the total charge on, say, the left edge is 
no longer a constant of motion: its time derivative defines 
the tunneling current It as follows: 



(r*[(o)^'fl,(o) - r*^'t^(o)*L(o)) . (6i) 



Looking back at Eq. H4()|l we see that the task at 
hand is that of calculating the correction to the displace- 
ment field propagator due to the interaction between the 
bosons. We will now show that this correction can be 
exactly expressed in terms of a tunneling current prop- 
agator and will provide a perturbative evaluation of the 
latter. 

First, let us introduce some compact notation. We 
define 



Bl = 



(62) 



where i = 1(2) for the upper (lower) component, and the 
associated phonon propagator 



-Q{t){Wn{t).B^^,]) 



(63) 
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Similarly we define 



(64) 



so that the phonon field propagator can be written as 

D{x^,x'p,t) = hi.^\,{x^)V^^,{t)^l{x'p) (65) 

(sum over repeated indices). 

The phonon operators satisfy the equation of motion 



where 



w„ 

—UJn 



In 



and 



(66) 



(67) 



(68) 



Then it is straightforward to verify that the phonon 
propagator satisfies the equation of motion 



where the auxiliary propagator 



(70) 



satisfies in turn the equation of motion 



with 



and 



{idt5,,-^^)gi(t) 



Mrit) ^ Mrit) - ^{HT)m 



-Mrit) (71) 



MT{t)^--Q{t){[lT{t),lT]) 

is the tunneling current propagator. 



(72) 



(73) 



This system of equations is readily solved by Fourier 
transformation with the following result 



X K 



(74) 



where [D^'^'^X^n' ^'^) noninteracting phonon prop- 

agator. Thus the tunneling correction to the phonon 
propagator is expressed in terms of the tunneling cur- 
rent propagator, as promised. 

We can now make use of this result to calculate the 
correction to the ideal conductances obtained in section 
III. Let us denote by cf^ the conductance obtained in 
the absence of tunneling and by 



5G, 



Gij " G. 



(0) 



(75) 



the correction due to the tunneling. After some straight- 
forward manipulations one arrives at 



SG,,^-^ lim V^„ai?(")](x„„0^;c^) 

(76) 

X [c.MT(c.)]^[i?(")]*(05,^-;^)e/3, 

where the indices 7 and S run over the two edges 
that are coupled by tunneling at a; = 0, and the 
Green's function of the noninteracting displacement field, 
[D<'°^]ap{x,x';uj), is given by Eq. 

As a concrete example, consider a four-terminal ge- 
ometry, as may be obtained from Fig. ^ by considering 
only terminals 1-4. Assume for simplicity that the mix- 
ing angle 9 is independent of x. Then from Eq. (|41|l we 
immediately get 



7 

iv 



(jj 

+ei-x) 



e(x) 



(77) 



-V 

u 



where the upper (lower) component refers to the left 
(right) edge and u — cosh0, v = sinh0. 



Substituting this in Eq. (|76l) we find 



(78) 
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where 



SGa(j{x, x') — — ie 



-20 



Van 



Mt{uj) 



iiJ~*i) LjJ 



uv u 

,,2 



e{x)e{-x') 

-Q{-x)Q{-x' 



-uv 

,2 



-uv V 

,,2 



e{~x)e{x') 

2 

uv 



V 

-uv u 



uv 

2 



(79) 



Putting this in Eq. (|78() and noting that = = 1, = ^,33 = ~1 (with the labels i, j as specified in the figure) 
we finally obtain the correction to the Landauer-Biittiker conductances of the ideal system: 



oGii ~ *e lim 



I uv 

,,2 



-uv —V 
0,2 



-uv —V uv V 



uv 

2 



(80) 



In Appendix D we show that 



_ Mt(^) 
gT = ^ lim = 4 



dt tlmG- (t) 



where G_ {t) = G_ (0, i; 0, 0) and 

G.{x,t-x',t') = ^: ^sl{x',t')^SR{x',t') :: ^-^^(0:, i)*i(x, t) 



(81) 



(82) 



Thus, the complete set of conductances has been expressed in terms of equilibrium averages of quasiparticle operators. 
Notice that the presence of a bias voltage Vq on the edge a modifies the time evolution of the corresponding 

quasiparticle operator from 'i!\,{xa,t) to '^l^{xa,t)e~'^ R^. The underlying physical assumption is, of course, that 
each edge is in equilibrium with a reservoir at potential V^. Under this assumption, the bias voltage dependence of 
the conductances can be calculated with no additional effort. 

To conclude this section we consider a specific experimental setup of Fig. 1."^^ The resistance Rxx of the quantum 
point contact is measured between terminals 3 and 4 in Fig. ^ 



V4-V3 
I 



(83) 



where / is the source-to-drain current. By considering that the constriction does not affect the source and drain 
probes, the full conductance matrix reads 



G., = 



/ 


1 














-1 


\ 







1 





-1 










-1 





1 + 5.911 


5gi2 


5gi3 


Sgii 












-1 + (5.921 


1 + <5.922 


6g23 


Sg2i 









-1 


'5.931 


(5.932 


1 + '5.933 


(5.934 




v 








'5.941 


<5.942 


-1 + '5.943 


1 + (5,944 


/ 



(84) 



As it is 



where the indices i, j run over {S", Z?, 1, 2, 3, 4} and the right bottom submatrix is given by Sgtj 
customary in the experimental setup we fix Vd — 0, Is = —Id = —I and h = I2 — I3 — Ia = 0. With these 
constraints, the equation (|29l) can be easily solved, and to the lowest non vanishing order in gT we get 



Rxx — 



e ^^(w -(- v)u gr 



-e cosh0 gT- 



(85) 



Notice that this perturbative result is valid only so long as Rxx is much smaller that the tunneling amplitude F 
must be sufficiently small for this to happen. 



VII. TUNNELING IN THE PRESENCE OF A CONSTRICTION 



Let us apply the formalism developed in the previous section to evaluate the resistance of a constriction of the type 
shown in Fig. The mixing angles 61 and 63 in the two external regions are assumed to be equal, while 62 {O2 > di) 
measures the strength of the interaction within the constriction. 
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The calculation of the correlation function G± can be perfomed by using the definition H57|l for '^^ and the Haussdorf 
lemma 



We start by considering the zero temperature limit where we obtain 



(86) 



G-(x, t] x' , t') — exp 



27re* 



A,fe>0 



(87) 



When we substitute the functions (p^'"'''\x) with those determined in eqs. 1)45(1 and 1(47(1 we obtain (we use x = x' with 
X inside the region of the constriction) 



G-[x^ t; x, t') = exp 



27re*' 



fe2>0 



where the coefficients , B", A'^ and B'^ are given by 1481) - ((501) . By assuming that the tunneling is localized only at 
the point x — x' = Q and substituting the expression 1(48(1 in this equation we obtain the key result 



G-{t) =exp 



47re*' cosh(26>2) ^ae^ 
Lve'^ cosh(26'i)'^ 



^ /cosh(26li2) -sinh(2e'i2)cos(fc2d)\ e'*'^'^^* 



fe2>0 



2sinr(e'i2)sin-'(fc2d) 



(89) 



where we have defined 9i2 ^ 0i — 02- For the function G+ (t) the calculation is similar and we can obtain G+ {t) from 
the above expression with the substitution t — > —t. 



Before going into the detailed analysis of the above ex- 
pression, it is useful to recall that in the limiting case 
012 = we recover the result of Wen^^ for the case of 
interacting edges 



G^{t) = exp 



47re*^ 
Lve^ 



-201 



fc>0 



iujkt 



k 



(90) 



Notice that the presence of the inter-edge interaction 
leads to a renormalization of the power-law behavior of 
the current-voltage characteristics via the factor e"^^^. 
The explicit form of the function G^ {t) can be obtained 
by using the well known analytical results 



Ecoa(nq) 1 , /„ „ / 
^^ = --ln(2-2cos(g)), 

ji=i 



n 



E 

n=l 



sin(nq) 



1 



(91) 



(TT-q). 



If g is a small quantity we have the approximate results 
cos(ng) ^ sin(ng) vr 



ln(g), 



E 

ri=l 



(92) 



In the case of the expression ((90(1 we can evaluate the se- 
ries, after defining the integer j as j = fciL/27r, obtaining 



G^(i) = exp 



2e* 



3-^"' In ( 



1-e^ 



(93) 



and in the limit of large system size ct/ L <^ 1 we have 



2Trict 

IT 



(94) 



where 5 assures the convergence of the series even when 
t ^ 0. This function is the propagator for the Luttinger 
Liquid model, with the anomalous exponent ^72-6"^^^. 
Notice that if we assume e* = ve we get for this ex- 
ponent 2ve~^^^ = 2j>. In this case the tunneling dif- 
ferential conductance at zero temperature is predicted 
to have a power law behavior with exponent given by 
2{e* /eY /v - 2. Let us go back to Eq.(|H51). First we no- 
tice that the additional fc-dependent factor in the sum of 
Ea. ((89() does not alter the logarithmic behavior at long 
times. To see this we define the quantity 



27r cosh(202) 
Tcosh(26'i) 

^ gik2C2t 



(95) 



^ k2 

k2>0 

cosh(26'i2) — sinh(20i2) cos(fc2d) 
1 + 2sinh^(6'i2)sin2(fc2d) 



and evaluate it numerically. To do this we calculate sep- 
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arately its real and imaginary part, 



(sin^(fc2C?)) =1/2 obtaining 



n o ,N cosh(26'i2) - smh(26'i2) cos(27rdn/i2) 

1 + 2 sinh^ 012 sin^{2TTdn/L2) 
cos{2TTnc2t / L) 



T r, 1^ i\ cosh(26'i2) - sinh(26'i2) cos(27rdn/L2) 
Imb-{t,d) — > 5 t: 



2sinh^ 6112 sm^2TTdn/ L2) 
sm[2nc2tn / L) 



(96) 



Notice that in these sums we have substituted k2 = 
21:11/ L2 where n is an integer and 



L2 = L 



cosh(26li) 
cosh(26'2) 



(97) 



takes into account the different speed of propagation of 
the waves in regions 1 and 2 (see Eq. H46|l ^. 

It is now useful to observe that the length of the con- 
striction introduces a characteristic time scale = d/c2, 
the travel time of an edge wave across the constriction. 
This clearly identifies a short (t < io) and long (t > io) 
time regime. In these two regimes, the equations H96|l 
may be approximated by taking the small and large d- 
limit in the A:-dependent factor. First we get for d — > 
the expressions 



ReS^{t,d^Q) =e 



4lnr2-2cos(^) 



2 

~-e-2ei2 In(27rc2</i2), 
ImS_{t,d^Q) =6-2912 (|(i_4c2t/L2)) 



In this limit the function G-(t) will then read 



(98) 



G-M^o{t) =[S- 



27ric2i 
L2 



(99) 



Remembering that the velocity C2 and the length L2 are 
related by Eq. 1)9 7|) we recover exactly the result one has 
when the constriction is not present. 

In the other limit d ^ 00 we have substituted in 
these two functions the averaged values, (cos(fc2d)) = 0, 



ReS-{t,d^ 00) =(1 + tanh2(6'i2)) 



1 , / „ „ /27rc2t 
— In 2 - 2 cos — 

2 V V ^2 
(1 +tanh2(6'i2))ln(27rc2t/L2), 



ImS-{t, d^oo) =(1 + tanh^(6'i2)) ( -(1 - 202^/^2) 



,{l + tanh^{9i2))-. 



(100) 



Again when we consider the function G_ (t) we get a 
power law of t 



2TliC2t 



''2(l+tanh''(ei2)) 



(101) 

In this case the presence of the constriction affects the 
exponent of this correlation function and can change the 
behavior of the tunneling amplitude. 

The two limiting regimes of short and long times of 
S-{t,d) are clearly visible in the full numerical evalua- 
tion as shown in Fig. 21 In the calculation of the real 
and imaginary parts of S^{t, d), we have fixed a value of 
d and then varied the value of C2t. As it is seen in the 
figures the two limits we have discussed are reached when 
d <C C2t 01 C2t. In Fig. ^we plot the numerical result 
for these functions as a function of C2t/ L for some value 
of d and fixed 612- We have restricted the calculation 
only to the limit of large system size ct/L, d/L ^ 1. The 
agreement of the calculated expressions with the approx- 
imate results (I98II . (jlOOf) is very good. Hence from now 
on we will use the simple expressions (|98|l and H100|) to 
carry out the calculation of the tunneling conductance. 

We then approximate the whole sum (EH) with a com- 
bination of two functions in the form 



S- (i, d) =e{t - to)S- {t, d^O) 

+ e{to-t){S-{t,d^ 



-A-) 



(102) 



where the two functions S-{t,d 0) and S-{t,d — s- 00) 
are determined by the corresponding limits for the func- 
tions ReS^ and ImS-. The factor A_ assures that 
S^{t,d) is a continuous function of t. With this ap- 
proximation we have separated the long time and short 
time behaviors of the response function. We then expect 
that the low energy behavior (which corresponds to the 
low bias voltage region) of the conductance will be dom- 
inated by the long time part of S-{t, d). Vice- versa, the 
response to a high bias voltage will be dominated by the 
short time behavior of S-{t^d). Within this approxima- 
tion the function G- reads 



G_(0 =e(t-io)G_,d^o(0 

+ e(io-t)G_,d^oo(t) 



X exp 



2e* 



-26, 



(103) 
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FIG. 4; a) Plot of ReS~{t,d) vs. In(c2t/-L2) for various val- 
ues of d. Observe the two different regimes for C2t > d and 
C2t < d. The two slopes agree very well with the approxi- 
mated result of Eg. 111)211 . We have chosen exp(Si2) = 1.5275 
in this calculation, b) Plot of ImS{t,d) vs. In(c2t/L2) for 
various values of d. We used the same parameters as a). The 
values for small and large C2t/ L2 agree well with the expected 
results (see Eas. (l98t and llOOt 'l. The downward curvature at 
large times arises from the finite size of the system used in 
the numerical calculation and disappears in the limit of large 
system size. 



Having obtained the expression for the function G- 
it is now possible to calculate the response function. We 
take into account the finite potential difference across the 
Hall bar via the replacement 



G-{t) G-{t)e' 



(104) 



where 



V1+V2 V3 + Vi 
Vt — 7: — Vh, 



(105) 



is the potential difference across the quantum point con- 
tact, and coincides with the Hall voltage. With this 
transformation, as is shown in the Appendix^ we get 



5t(wt) = 4^^ ^Im / dte'^^^HmG-it). (106) 

h-^ duJT Jo 



A lengthy but straightforward calculation gives the ex- 
pression for the tunneling conductance at zero tempera- 
ture (see Appendix IeI for details) 



5t(wt) 



4|r|2e*^to\ / a 



ft3 



sm 



/"Ka\ 



duJT 



TtoT ^ (cos sign(wTio)Rer(l - a, -iurta) + sin (Jy') ^^^^^(l - a, -iujTto) 



cos 



^ ) sign(wTio)(r(l - /?) - Rer(l - /?, -nurto)) - sin ( ^ ) Imr(l - /3, -lajrio) ] 



(107) 



where we have defined 



n 

a — 

1^ 



-261 



(108) 



2 e* 



^e-'''ni + tanh"(0i2)). 



a is a short-distance cut-off, and r(zi,Z2) is the incom- 
plete r function^. 

In Fig. Owe plot Rxxi^^r) in the case that the inter- 
edge interaction is confined to the region of the constric- 
tion (i.e., we set 6'i = ^3 = and let 62 assume several 
different values) . Experimentally, 62 can be increased by 



narrowing the constriction by the application of a gate 
potential. When 92—0 there is no interaction and Rxx 
diverges as V^~^ at low bias. This low-bias behavior does 
not change upon increasing 62 because the long time be- 
havior is dominated by the exponent a which does not 
depend on 6*2. At larger bias voltage on the other hand. 



Rxx behaves as Kj 



(3-2 



Furthermore, the plot of Rxx 



shows oscillations, which become more pronounced with 
increasing 62- We can express the period of these oscil- 
lations in terms of the physical parameters of the theory 



AVt = 



h 



hci cosh 26*1 



(109) 



e*to e*d cosh 26*2 
The frequency of the oscillations increases with increasing 
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FIG. 5: Plot of the resistanc e Rx x/\V\^ given by Eq. ||H3 
with the qt calculated in Eq. 110711 for various values of 62 at 
fixed 6\ = 0. The oscillations at large bias voltage becomes 
more and more pronounced with increasing 62- 



02 as it is apparent in Fig. |S| 

The finite temperature behavior of the tunneling re- 
sistance can be derived from the zero-temperature be- 
havior of the same quantity by means of the conformal 
transformation— 

Notice that we are using units in which h — ks = ^, 
where fcs is the Boltzmann constant. The correct phys- 
ical dimensions are restored via the substitution T — > 
ksT/ h and this is understood in the Eqs. (|lll|l and H112|) 
below. Making the transformation (|110|) in Eqs. H99|) and 
(|101|l . and substituting in Eqs. (|102|l and (|103|l we obtain, 
after lengthy calculations (see the Appendix IeI for more 
details), 



5r(t^T,r) =4- 



irr 



\2TrcJ 



aT 
"2^ 



sin(^) d ^ 
— ^ ^ ' Im 

sinh"(7rrto) dujx 



LOT 



1-/3 



a ~ I 



■UJt_ 

T 



F a,l;l 



a 
2 



LOT 



r 



/3,l;l-,- 



LUt 



2ttT' 1 



(111) 



where F is the hypergeometric function of four arguments 
(also indicated as 2F1) and B the Euler beta function—. 
In the case Oi — 62 we have a = f3, the first and third 
term cancel against each other and we recover Wen's re- 
sult 



5t(wt) 



iri 



d 

duji 



( a \ faT 
\2ttc) \ c 



-ImB 



'a 
.2 



I , 1 

2t:T 



(112) 



In Fig. Efa-d) we plot the differential resistance Rxx vs. 
bias voltage for a system without inter-edge interaction 
(dashed line - 6*2 = = 0) and with inter-edge interac- 
tion (solid line - 0i = 0, ^2 = 1) for different values of 
irdT/ci = 0.1,0.5,1, and 1.5. The non vanishing value 
of 02 within the constriction induces oscillations in the 
Rxx vs. LOT relation with the same period as in the zero 
temperature case. However, we now have a maximum at 
zero bias voltage and two minima at finite bias voltage. 
This behavior is due to the fact that the temperature in- 
troduces a new energy scale. When the e*V > ksT we 
are essentially in the zero temperature case and the re- 
sistance Rxx decreases with decreasing bias voltage (see 
Fig. O. But, when the e*V < ksT the resistance turns 
around and begins to increase, reaching a maximum at 
zero bias. This behavior implies the presence of two min- 
ima located at bias voltages of the order of magnitude of 
ksT/e*: these are clearly seen in Fig. El The finite 
value of Rxx at zero bias (independent of Vt to first or- 




FIG. 6; Plot of the differential resistance -Ra:i/|rp vs. lot for 
a system with inter-edge interaction within the constriction 
(continuous line, 61 — 0, O2 — 1) and without inter-edge inter- 
action (dashed line, 6*1 = ^2 = 0). The four curves correspond 
to different temperatures: vrTd/ci = 0.1 (a), 0.5 (b), 1 (c), 
and 1.5 (d). 



der) indicates that the constriction is behaving like an 
ohmic resistor in this regime, even though the resistance 
is strongly temperature-dependent. 

The presence of a constriction adds another energy 
scale in the problem, associated with the inverse of the 
characteristic time to- For temperatures smaller than 
h/to the low bias behavior is dominated by the same ex- 
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ponent a (cf. Fig. |B| (a,b)) irrespective of whether the 
inter-edge interaction is present or not. When the tem- 
perature, instead, is greater than h/t^ the exponent /?, 
which depends on the strength of the interaction within 
the constriction, controls the behavior of R^x (cf. Fig. 
Elcjd)). As a consequence the minima at finite bias are 
generally deeper and shift to lower voltages. 

The effect of the constriction depends quantitatively 
on both the inter-edge interaction parameter O2 and the 
temperature. To appreciate this we plot in Fig. |3the dif- 
ferential resistance Rxx for different values of the inter- 
edge interaction and the temperature. More specifically 
we have plotted Rxx without interactions (^2 = = 0) 
and with interactions within the constriction (^i = 0, 
6*2 = 0.2) for TidT/ci = 0.5, 1, 5, and 10. We notice that 



range. 

At lower temperatures, on the other hand, the experi- 
ment shows a completely different behavior which is not 
qualitatively consistent with the present theory irrespec- 
tive of the presence of inter-edge interactions. Strong 
tunneling effects^, which can be treated by the ther- 
modynamic Bethe Ansatz, are not likely to explain the 
unexpected decrease in Rxx that is seen at these tem- 
peratures. This clearly suggests that a different physical 
mechanism comes into play at these temperatures and 
some additional physical input is needed. One could, for 
instance, speculate that, within the constriction, the hy- 
drodynamic approximation may be too crude and better 
treatment of the edge structure may be required. This 
is, however, outside the scope of the present work. 
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FIG. 7: Plot of the differential resistance i?n/|rp vs. fre- 
quency with and without inter-edge interaction within the 
constriction. Solid line - 0i = 0, = 0.2; Dashed line - 
9i = 62 ~ 0. Temperatures are ndT/ci — 0.5 (a), 1 (b), 5 (c), 
and 10 (d). 

the effect of the inter-edge interaction disappears at suffi- 
ciently low temperature, since it is always the long times 
exponent a that matters in that regime. The effect of 
the interaction shows up upon increasing the tempera- 
ture above the crossover energy h/to: the latter decreases 
with increasing 62- Such a trend is clearly seen by com- 
paring Figs, inland [71 We note that similar crossover ef- 
fects in the temperature and voltage behavior have been 
discussed also in the context of transport in quantum 

Finally we would like to comment about recent 
measurements of tunneling characteristics through a 
constriction^ in the weak inter-edge tunneling regime 
at high magnetic field. At relatively high temperatures 
(T > 400mK) the experiment clearly shows the emer- 
gence of a zero bias peak in the differential longitudinal 
resistance which is qualitatively consistent with the re- 
sults presented above. The experiment also shows well 
defined minima at finite bias voltage, which, according 
to the previous discussions may reveal the effect of the 
constriction. In fact, the system without inter-edge in- 
teractions never shows deep minima in this temperature 



VIII. CONCLUSION 

In this paper we have extended the derivation of the 
^LL model to arbitrary values of the filling factor v. 
We have developed a theory to calculate the Landauer- 
Biittiker conductances for various experimental setups, 
taking into account both inter-edge tunneling and inter- 
edge interactions. In the absence of tunneling, our model 
recovers the usual fractional Hall conductance, even when 
an inter-edge interaction is present. The breaking of 
translational invariance, due to the constriction, does not 
change the low-frequency behavior of the conductance as 
long as tunneling can be neglected. 

We have then discussed the effect of inter-edge tun- 
neling. Tunneling destroys the exact quantization of the 
Hall conductance. We have calculated the tunneling con- 
ductance (related to the resistance of the constriction) to 
the second order in the tunneling amplitude. A problem 
with the present form of our theory is a fundamental un- 
certainty about the value of the effective charge e* of the 
quasiparticles at generic filling factor v. This remains a 
major open theoretical question. 

The presence of the constriction introduces a finite 
time-scale (the time it takes an edge wave to travel along 
the constriction) and gives rise to different short- and 
long-time behaviors of the tunneling propagator. The 
long time (low frequency) behavior is dominated by an 
exponent that coincides with the one well known in the 
literature. The short-time behavior is dominated by 
a different exponent, smaller than the long-time expo- 
nent. The interplay between the two exponents intro- 
duces small oscillations in the tunneling conductance 
which can possibly be used to measure the amplitude of 
the inter edge interaction and the velocity of the modes 
in the constriction. 
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APPENDIX A: COMMUTATION RELATIONS 

To derive the commutation rule of Ea. (|ll|) . one integrates Eq.® with respect to y and y' according to the pre- 
scriptions given in Eq. H10() . The delta function makes the commutator non vanishing only for points belonging to the 
same edge and the final result will have a factor Sap- We get 



[8p{xa),5p{x'a)\ = i^^y AyAy'{dx^po{xa,y))dy5{xa- x'a)5{y - y') 

-U^ dydy'{dyPo{xa,y))dx^S{xa ~ x'a)S{y - y'). (Al) 



One observes that the first term on the right-hand side vanishes, while the second term gives, after making the 
integration over y' and y, for a ^ L 

[6p(xL),dp{x'^)] = ~ie'^(po{xL,d) - pn{xL,0))d.x^S{xL - x'^) 

= -ifpoix)d,JixL-XL) (A2) 

and for a = i? 

[Sp{xR),Sp{x'jj)] = -if{po{xR,0)-po{xR,-d))dxaS{xR-x'j^} 

= ie^Po{x)dxJ{xR-x'ji) (A3) 

where d indicates the distance from the edge at which the density has reached its bulk value, pa{x). The different 
order of the limits of integration for the two edges gives the relative minus sign between the edges. 



APPENDIX B: PROPERTIES OF THE 
EIGENVALUE PROBLEM 

In this appendix we want to study some analytical 
properties of the equation l|18() . First of all let us de- 
fine the operators 

Ma = isad^^, (Bl) 
I' f°° 

Ha,f3 = dx'p dx^V{xa,x'p)dx,^. (B2) 

With this definition we rewrite the equation of motion 
(|18|l in the compact form 

ujMip^Hip. (B3) 

It is easy to see that H and M are hermitian operators 
and we request that H is positive definite (this will assure 
the stability of the physical system). 
Let us define the auxiliary function 

^ = H^ip (B4) 
which is a solution of the equation 

= (^H-^AIR-^^^ ^ Am (B5) 

if (/J is a solution of H18|) . Because M is a hermitian 
operator we have the results: 



1. the set {^'} of solutions forms a complete base of 
the Hilbert space, 

2. the orthonormality condition is 

/ dXaKiXa)'i>miXa) = Sn,m, (B6) 

3. the completeness relation 

*„(a;a)**„(a::^) = Sa^pSix^ - x'p). (B7) 

n 

Because there is a one-to-one relation between and 
^ we have the following properties of the solutions of 
equation lfTH|) : 

1 . the solutions Lp form a complete base of the Hilbert 
space, 

2. they are orthogonal with respect to the scalar prod- 
uct 

(</?n,</'m) = X! / UJ.aiPn{,Xa)MaPm{.Xa), (B8) 

3. they satisfy the completeness relation 

-iYu!nfn{Xa)'Pn{x'f3)Si3dx'^ = Sa,l3S{Xa - x'p) . (B9) 
n 
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We obtain the relations reporte d in t he text if we nor- 
mahze the functions as ipnl \f\^\- 

Now we want discuss the degeneracy of the eigenvalues 
of the equation H18(l : 

• If ^Pm{?^a) is a solution with given eigenvalue 
then the function (^*j(xq,) is also a solution with 
eigenvalue ct;_„ = — 

• If ipm,{xa) is a solution with given eigenvalue ajm 
then the function f3^"i{xi3) is also a solution 
with eigenvalue —iVm- 

Then we have that if ipm{xa) is a solution then 
a^l^ip*{x[3) is still a solution with the same eigenvalue: 
that is the solutions of problem l|18(l are doubly degener- 
ate. 



APPENDIX C: CONSERVATION LAWS 

In the case Va,f3{x~x') = Va,i3{x)S{x — x') the quantity 

fiix)^^,l3Vl3ix) = (p[{x)(pL{y) ~ ip''jiix)(pR{x) (CI) 

is conserved 

d. Hix)a^.j3M=^)) = 0- (C2) 

The proof of the existence of this conserved quantity rests 
on the basis of the existence of the inverse of the matrix 
Va^[3{x) for every value of x. Consider the equation of 
motion and its complex conjugate for the displacement 
field wave function (p{x) 

iuj(fa{x) ^ -^Vafi{x)d^tpp{x) (C3) 
ZTT 



-iujipl,{x) 



VSa 

27r 



d,^}{x)Vp^{x). (C4) 



The conservation law follows by first taking Vap on 
the left-hand side of both equations and then multi- 
plying the first (second) equation on the left (right) by 
^ff'^^ai'^afj'Pp)^ finally summing the two equations. 



APPENDIX D: CALCULATION OF THE 
TUNNELING PROPAGATOR 

In this Appendix we derive Eq. 181|) . The first task is to 
compute MT{t). We do it to second order in perturbation 
theory in F. The first term in the definition of Mxit) 
propagator. Since it is already second order in F, we only 
need to evaluate its average in the unperturbed ground 
state. This can be expressed in terms of the correlation 
functions G^{t';t) defined in (|82|) obtaining 

([HT{t'),HT{t)]) = Ai\T\HmG.{t';t) (Dl) 



Notice that in this expression we have dropped the 
anomalous averages that appear when one considers the 
average value of several field operators This is jus- 
tified by the presence of the fermion operator U in the 
definition of the quasi-particle operators (|57|l . Hence the 
contribution of tunneling current propagator to Alxit) 
reads 



:/T(t),/T(o)] 



-Ai—\T\HmG-{t;0). 



(D2) 



We now consider the other term in Alxit). This is the 
average of the tunneling Hamiltonian and is only first 
order in F so that we need to compute the first order 
correction to the ground state as well. We get 



dt' {[HT{t'),HT{t)]) (D3) 
dt' 4i|FpImG_(t';t) (D4) 



We are now ready to compute the Fourier transform of 
MT{t). We get 



lim 



4e*2|F 



lim 



i lim 



dt 





cos(a;t) 



UJ 



-ImG- (t) 



.4e 



*2|p|2 



ft3 



I dttlmG^it). 
Jo 



(D5) 



The presence of a voltage difference between the edges 
can be taken in to account by means of the transforma- 
tion 



G-it) 
G+{t) 



iujTt 



G-{t), 
G+{t). 



(D6) 



from which we get 

'[HT{t'),HTit)] 

:/T(t),/T(0)]\ = -4^ 



4i|Fpcos(wTt)IniG_(t';t), 
^*2|-ri2 

■ cos(cjtOI™G- (*) 



(D7) 



ft2 



and we finally arrive to the expression for gT 



i lim 

uj—>0 



Mt{uj) 



4e*lF|- 



UJ 



fi3 



duJT Jo 



dt e*"^*ImG_(i). 

(D8) 



APPENDIX E: EVALUATION OF INTEGRALS 

In this appendix we provide a few details concerning 
the evaluation of the integral occurring in the calculation 
of the Fourier transform of the response function G_ (t) . 
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t„ 



X 



FIG. 8: The integration path used to evaluate the integral 
llEll l when the frequency is positive. A similar path, closed 
in the lower plane, is used in the case uj < 0. 



In the zero temperature case, we must evaluate an inte- 
gral of the form 



dt{S±it) 



(El) 



where a is a positive real number. To do this we go in the 
complex plane of the variable t and consider, for positive 
frequency uj, an integration path as the one shown in Fig. 
IHl A specular path in the lower half plane must be used 
for negative frequency. We observe that the integrand 
function has no poles in the complex half-plane of t with 
a positive real part. Hence the integral on the whole path 



is zero. The integral on the arc vanishes by letting the 
radius go to infinity. 
As a result we get 



dt{5 ± itye'"^* i(Tl)~"w""^r(l - a, -iuto). 

(E2) 

The case = can be carried out by calculating the 
convolution product between the Fourier transform of the 
0(i) function and the integral 



/oo 
dt{S ± itye'^' =2 sin(7ra)e*^ " l^^l^i-" 
-OO 

xr(l + a)(Tl)-^^", 



(E3) 



obtained by cutting the complex plane along the imagi- 
nary axis, starting from t — ±i(5. 

In the finite temperature case, we need to calculate the 
integral 



dt e" 



sinh(7rTi)" ' 



(E4) 



One can easily obtain the result reported in the text by 
means of the substitution s — e"^'^-^* which reduces the 
above integral to the definition of the hypergeometric F 
function of four arguments^SiS. The integral with to = Q 
can be easily obtained in the limit ^ by using the 
corresponding limiting expression of the hypergeometric 
function F in terms of the Euler beta function. 
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